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This paper discusses the optimum design of discrete, directive antenna 
arrays of arbitrary geometrical configuration in space, when the excitations 
and spatial positions of the elements vary in a random fashion about their 
nominal values. Under certain assumptions the expected power pattern of an 
array turns out to be the power pattern of the nominal array, plus a "back- 
ground" power level which has the same dependence on direction as the pat- 
tern of a single element. A set of excitations which maximizes the theoretical 
directivity of an array may correspond to a superdirective design, in which 
the background power level will completely swamp the desired pattern unless 
the excitations and positions of the elements are controlled with extraordi- 
nary precision. A method is given for maximizing the gain of the array while 
holding the expected background power level constant, when the precision 
with which the excitations and positions can be controlled is known. The 
method is illustrated with numerical examples. 

1. INTRODUCTION 

The effects of random variations on antenna patterns have recently- 
been discussed in a number of papers which treat more or less special 
cases. Several authors 1, 2 ' 3 are concerned with linear Chebyshev-designed 
broadside arrays, and in particular with the effects of manufacturing 
variations on the patterns of slotted waveguide arrays. Ruze has de- 
rived the expected pattern of a plane, rectangular array in which the 
positions of the radiators are rigidly fixed, while only their amplitudes 
and phases are variable. Less attention seems to have been given to the 
possibility that the positions of the elements may also vary in a random 
fashion, even though this latter situation can arise in any electromag- 
netic or acoustic array whose elements are not rigidly supported. 

1 L. L. Bailin and M. J. Ehrlich, I.R.E. Trans., PGAP-1, pp. 85-100, Feb., 1952. 
» D. Ashmead, I.R.E. Trans., PGAP-4, pp. 81-92, Dec, 1952. 

3 H. F. O'Neill and L. L. Bailin, I.R.E. Trans., P(JAP-4, pp. 93-102, Dec, 1952. 

4 J. Ruze, Nuovo Cimento, 9, Supp. 3, pp. 304-380, 1952. 
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This paper derives a relatively simple statistical expression for the 
expected power pattern of an array when the excitations and positions 
of its elements are subject to independent random variations. The re- 
sult is not restricted to linear or even to plane arrays, but is valid for 
arrays of arbitrary geometrical configuration in space. In the special 
case when the excitations are variable and the positions of the elements 
are either fixed or subject to random displacements with a spherically 
symmetric distribution, the expected power pattern is the pattern of the 
nominal array, plus a "background" power level which has the same 
dependence on direction as the pattern of a single element. The expected 
background power level is proportional, for small errors, to the sum of 
the mean-square errors in excitation and position. 

We consider in particular the application of our results to the prob- 
lem of designing superdirective arrays. A superdirective array is one 
having a beamwidth in radians much less than the reciprocal of the 
largest dimension of the array in wavelengths. It is well known that such 
narrow beams can be designed on paper, but only by employing heavy 
cancellation between adjacent elements. If the excitations and positions 
of the elements of a superdirective array are not controlled with great 
accuracy, the background power level due to random errors will com- 
pletely swamp the desired pattern. This corresponds to the familiar fact 
that a really superdirective array has to be constructed with extraordin- 
arily high precision in order to give anything like its calculated perform- 
ance. 

A method is given for computing the excitations which maximize the 
gain of an array of any specified geometrical configuration, while hold- 
ing the expected background power level due to random variations 
constant. It is assumed that the excitation coefficients can all be con- 
trolled to the same per cent accuracy, and that the element displace- 
ments, if any, are distributed with spherical symmetry. If the random 
variations are taken to be zero, or if no restrictions are placed on the 
background power level, the present procedure becomes equivalent to 
the methods which have been described in recent papers 5, 6 on the max- 
imum gain of an arbitrary array. 

An interesting result of the analysis relates to arrays in which the 
elements are all excited with equal amplitudes, and with phases such 
that their fields add in phase in a specified direction in space. This ar- 
rangement may justifiably be called the normal excitation, since it is 

5 A. I. Uzkov, Dokl. Akad. Nauk SSSR, 53, pp. 35-38, 1946. 
« A. Bloch, R. G. Medhurst, and S. D. Pool, Proc. Inst. Elect. Engis., 100, Pt. 
Ill, pp. 303-314, 1953. 
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often adopted in practice as a means of obtaining a beam in the desired 
direction. It is proved that, of all possible excitations of the array which 
produce a main lobe in the given direction, the normal excitation leads 
to a pattern which is most insensitive to the effects of random errors. 
Furthermore, if all the elements have the same ohmic resistance, the 
normal excitation produces the highest power flow in the direction of 
the main beam for a given rate of heating the elements. 

It is not invariably true that maximum gain is incompatible with 
minimum sensitivity of the pattern to random errors. The normal ex- 
citation maximizes the gain of any array of isotropic elements in which 
the distance from every element to every other element is an integral 
number of half wavelengths. However, for most arrays the gain can be 
increased above that obtainable with the normal excitation, at the ex- 
pense of a (possibly enormous) increase in the sensitivity of the pattern 
to random errors; and some sort of compromise will have to be struck. 

Several types of symmetry are commonly found in antenna arrays, 
and some of these symmetries force very simple relationships to hold 
among certain of the optimum excitation coefficients. We discuss the use 
of symmetry to reduce the amount of computation necessary in design- 
ing an optimum array. As an example, these considerations are applied 
to the design of an array of four elements located at the corners of a 
tetrahedron, and also to a four-element end-fire array. Curves are ob- 
tained which illustrate the relationship of gain to pattern sensitivity for 
several such arrays of different dimensions in wavelengths. 

It is shown that if an arbitrary array of isotropic elements is excited 
successively to have maximum directivity in different directions, the 
average value of the maximum directivity over all directions in space is 
equal to the number of elements in the array. The excitation required 
to produce maximum directivity will naturally depend on direction, and 
considerations of pattern sensitivity are ignored. The significance of 
this result is that if an array configuration permits an abnormally high 
gain in a certain direction, as is theoretically possible, for example, with 
a very short end-fire array, then there exist other directions in which 
the maximum gain of the same array is abnormally low. 

2. STATISTICAL FORMULATION 

Consider an antenna array of n elements. We shall call the elements 
radiators, though they may equally well be thought of as receivers. We 
assume that each element has the same directivity pattern s(u) with 
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respect to a fixed set of axes,* where u is a unit vector representing a di- 
rection in space, and s(u) is a complex-valued vector function giving the 
amplitude, phase, and polarization of the radiation field over a large 
sphere centered at the element. For acoustic fields, s(u) is a scalar 
function. The average density of power flow in the direction u is propor- 
tional to 

s 2 (u) = s(u)-s*(u) (1) 

where an asterisk denotes the complex conjugate quantity. 

Let the excitation of the fcth element be A k . The complex numbers 
Ai A 2 , • • • , A n will be called the excitation coefficients of the array. 
Let the position vector of the fcth element relative to an arbitrary origin 
be R fc . The field strength produced by the array at the point at the end 
of the vector Ru from the origin will, for large R, be proportional to 

s(u)/(u)/i2 (2) 

and the average density of power flow at the same point will be propor- 
tional to 

«*>(u)/fl 2 (3) 

where the power directivity pattern $(u) is 

*(u) = s 2 (u)|/(u)| 2 (4) 

and the array factor /(u) is 

/(u) = £ A k exp (#R*-u) (5) 

As usual, /3 denotes 2t divided by the wavelength X. 

Let us assume now that the excitation coefficients and the positions 
of the elements actually have some random scatter about their mean or 
expected values. We shall calculate the expected values of the field and 
power patterns. These expected values may be regarded as averages 
taken over a large number of different arrays, or they may be thought of 
as long-term time averages for a single array whose parameters vaiy 
with time in a random fashion. We can adopt the latter point of view 
when dealing with an array whose elements are not rigidly intercon- 

* No difficulties would result in the statistical analysis from assuming a differ- 
ent pattern for each element, hut the added generality would complicate some of 
the following work, and it is unnecessary for the great majority of practical arrays. 
We also ignore the possibility that the orientations of the elements might be 
subject to random variations. 
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nected, but are subject to displacements relative to one another in the 
course of time. 

To make matters precise, assume that the excitation coefficients are 
given by 

A k = a k + a k (6) 

where a k is the expected value of A k and the a k s are independent random 
complex variables with mean zero. Also let 

R* = r* + Q k . (7) 

where r* is the expected value of the position vector R k and the p*'s 
are independent random vectors with mean (0, 0, 0), all having the same 
statistical distribution. 

We can now write down the expected values of the field and power 
patterns. Denoting expected values by angular brackets, we find for the 
field strength, 

(s(u)/(u)> = s(u) £ (A fc )(exp (ipRk-u)) 

" (8) 

= s(u) X) a k exp (ipr k - u)(exp (i/3p*-u)> 

= (exp (?/3p-u)) s(u)/ (u) 

where q is a random vector having the same distribution as the Q k 's, 
and /o(u) is the nominal array factor 

/o(u) = X Ok exp (f/3r fc -u) (9) 

which results when the excitation coefficients and positions all have 
their expected values. 

The norm of the general array factor may be written 

n n 

I /(«) I 2 = Z A k exp (ipRk-u) 2 Aj* exp (- t/3R r u) 

= 53 YJ fak + a k ) (aj* + aj*) exp [^(r* - ry) • u] exp [ip(g k - qj) ■ u] (10) 

+ E(Ot+ a*)(0»*+ at*) 

where the primed summation sign indicates that the term for which 
j = k is to be omitted. Taking expected values and recalling that the 
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random variables are independent, we obtain 

<l/(u)| 2 ) 

= Z £' o-kdj* exp [ip(r k - tj) -u] (exp (i/3p*-u) > (exp (-ij8p,-u) > 
fc-i 1=1 

+ Eki 1 + i;<i«*i 2 ) 

fc=l *=1 



(11) 



= I (exp (%p 9 -u) ) | 2 | /o(u) | 2 + E <| a, | 2 > 

+ [1 - | <exp(^p-u)>| 2 ]i:|a,-| 2 

where the last step follows by adding and subtracting the terms with 
j = k which were omitted from the double sum. 

Multiplying through by the power pattern s 2 (u) of a single element 
gives the expression for the expected power pattern of the array, namely 

n 

<$(u)> = I (exp (#p-u)> | 2 *»(u) + « 2 (u) E (I «* T> 

(12) 

+ [1 - | (exp (tfo-u) > |V(u) Z I o fc I 2 



where the power pattern of the nominal array is 

$o(u) = s 2 (u)|/o(u)| 2 (13) 

3. SPECIAL CASES 

If the positions of the elements of the array are supposed to be exactly 
known and rigidly fixed, then the displacement vectors g k are identically 
zero, and the general result derived above reduces to 

<*(u)> - *o(u) + s 2 (u) E (| a, | 2 ) (14) 

k=l 

Equation (14) has a simple physical interpretation. It asserts that 
the expected power pattern is the power pattern of the nominal array, 
plus a "background" power level which has the same dependence on 
direction as the pattern of an individual radiator, and is proportional to 
the sum of the mean-square errors of the excitation coefficients. Of course 
for any particular array the background power will not have exactly the 
directional dependence s 2 (u), but will exhibit fluctuations depending on 
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the particular set of errors in the excitations. 7 However, in order to have 
the over-all pattern be a good approximation to the nominal pattern 
<Pn(u), it is necessary to hold the expected value of the background 
power well below the maximum value of «t>o(u). We shall discuss the im- 
plications of this requirement presently. 

If the displacements of the elements are not identically zero, then 
we denote the cartesian components of the vector g by £, 77, f , and their 
joint probability distribution by P(£, 77, f ). P(£, 77, f ) df d-q f/f represents 
the probability that the end of the vector g, drawn from the origin, will 
lie in the volume element r/£ dr\ d$ centered at (£, 77, f). Then the ex- 
pected value of exp (i0g-u) is given by 

(exp (tj8p-u)) 

= / / / exp [#({ti, + yuy + fuj] P(€, V, f) di dv <fr 

where w x , •»/„ , u z are the components of the unit vector u. 

In order to evaluate the integral (15) for any particular array, one 
should choose the most plausible joint distribution function P(£, 77, f) 
that his physical insight permits. In general (exp (i/3p-u)) will depend 
upon the direction of u; but if the distribution of {> is spherically sym- 
metric,* the expected value will be independent of direction. In this case 
P(.£) V, f) is a function only of the magnitude p of the displacement 
vector. To evaluate the integral (15), take the f-axis parallel to u, and 
let p(p) dp be the probability that the length of the displacement vector 
lies between p and p + dp. Then 

(exp (i/3p-u)) — — / / / exp (ifip cos 6)p(p) sin ddpdd d<p 

47T Jo JO Jo 

(16) 
= [ [(sin j8p)/j8p] p(p) dp = <(sin /3p)//3p) 

Jo 

If the cartesian components of g are assumed to be normally and 
independently distributed with mean zero and equal variances «r 2 /3, 
so that the root -mean-square value of p is a, one has 

p(p) = 3(6/tt)"' 2 (p 2 A 3 ) exp (-3p72<r 2 ) (17) 
and 
(exp (iffg-u)) = exp (-&V/Q) (18) 

7 The statistical distribution of those fluctuations has been discussed by Ruzc 
in Reference 4. 

* Under some circumstances the displacement vectors may be constrained to 
lie in a plane. If their distribution is circularly symmetric, and if we confine our 
attention to directions u lying in the plane of the displacements, the results will 
be formally similar to those obtained for the spherically symmetric case. 
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For any spherically symmetric distribution of g Ave may define a 
parameter 5 2 by 

8 s = | <exp0/3e-u)>r ~ l ( 19) 

If <r, the root-mean-square value of p, is small compared to the wave- 
length, then 

o 2 « /3V/3 (20) 

for the normal distribution (17), and also for other distributions which 
taper off for large p with comparable or greater rapidity. 

From equations (12) and (19) we obtain the normalized expected 
power pattern for a spherically symmetric distribution of displacements, 
namely 

u+y) <*(«)> 

r "1 ( 21 ) 

= * (ii) + s 2 (u) (1 + 5 2 ) E (| ot k | 2 > + 5 2 E I °>k I 2 

Again the expected pattern turns out to be the nominal pattern plus a 
background level with the same distribution as the pattern of a single 

element. 

In what follows we shall idealize the problem somewhat by assuming 
that the excitation coefficients A k can all be controlled to the same rela- 
tive accuracy. 8 Precisely we suppose there is a small number e such that 

<K| 2 )= e s | a* I s , fc = l,2, ••• ,n (22) 

Thus (21) becomes 

(1 + 8 2 ) (*(u) ) = *o(u) + AV(u) g | a k | 2 (23) 

where 

A 2 = (1 + 5 2 )e 2 + 5 2 W 5 2 + * (24) 

and the last approximation is valid if 5 2 is small compared to unity. 



8 Ruze's assumptions in Reference 4 amount to taking | A k \ and arg A k as 
independent random variables with means I k , d h , where the coefficients^ of the 
desii 
tributed 




proportioi.c. »„ , -* , . 

looks more complicated than equation (14) mainly because the expected value ot 
A k differs from h exp (idt) by a constant factor. 
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4. CRITERIA FOR GOOD PATTERNS 

The preceding statistical analysis is equally valid whether the nominal 
pattern $o(u) is to have some specified beam shape for a particular appli- 
cation, or merely to provide as narrow a beam as possible. However the 
latter case is of greater practical interest, and so we shall consider hence- 
forth only the design of highly directive or pencil beam arrays. 

We suppose that the number and configuration of the radiators are 
fixed in advance, and we have to choose the excitation coefficients to ob- 
tain the desired pattern of the nominal array. We shall denote the set 
of n complex numbers (a } , a 2 , • • • , a„) by the single symbol a. Our 
problem is to find a way of exciting the array so as to produce a narrow 
beam in some direction Uo , i.e., we must choose a in such a way as to make 
<p (1 (uo) small for all directions u not near u . When we try to formulate 
the problem more precisely than this we find a large number of alterna- 
tives. 

As a mathematically tractable criterion for a good pattern, we shall 
stipulate that a is to be chosen so as to minimize the generalized gain 
function 

, . 3>o(Uo) 

G{a) = T~7 (25) 

4tt 



~ f <p (uV(u) da 



where iv(u) is a non-negative weight function which may be chosen at 
pleasure, rffi is an element of solid angle, and S is the surface of the unit 
sphere. Since «$>o(u) is proportional to the density of power flow in the 
direction u, G(a) represents the ratio of power flow in the direction u n 
to a weighted average of power flow in all directions over a sphere. If 
one were interested in the radiation pattern only in some plane contain- 
ing Uo , the weighted average could be taken in all directions around a 
circle with only minor changes in the formal analysis. 

A few comments may clarify the significance of the function w(u). If 
itf(u) is taken to be identically unity, then G(a) is just the conventional 
gain of an arra}' with the set of excitations a. However, it is well known 
that merely maximizing the gain does not alwa} r s produce a pattern 
with low side lobes. If it is important to prevent large fields from being 
radiated in specified directions, one may choose u>(u) to be unity over 
the set of unwanted directions and zero elsewhere; and in principle even 
more complicated choices of w(u) could be made to discourage side 
lobes. For receiving arrays, tv(u) may similarly be chosen to discriminate 
against the reception of spurious signals from particular directions. 



G46 
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In terms of the excitation coefficients, the function G(a) has a simple 
formal expression, namely 

2 

2 at exp (i/3r*-u ) 



G(a) = s 2 (u ) 



(26) 



2 Zl hj k a k a,j' 

k=l ;~1 



where 



fc» = _L f exp [#(r* - r y )-u]s 2 (u)w(u) dS2 (27) 

47T Js 



and obviously 



h kj = hj k * 



(28) 



The coefficients h jk depend only on the weight function, the pattern of a 
single element, and the positions of the elements (in terms of wave- 
length). When the elements are isotropic and the weight function is iden- 
tically unity, we have 

sin fifjk 



h ik = 



0r y * 



(29) 



where r jk is the distance from radiator j to radiator k. In any case 0(a) 
is the quotient of two Hermitian forms in the excitation coefficients, and 
can therefore be maximized by standard mathematical techniques. 

It may be noted that the use of an integrated criterion such as (25) 
for keeping the fields small away from the main beam does not absolutely 
guarantee the absence of undesirably high side lobes in particular direc- 
tions for any given array. To be sure of keeping all side lobes below a 
certain level, we should choose a to maximize some such expression as 

*° (U0) (30) 



T(a) = 



max 3>o(u) 

Q 



where the denominator is the maximum value of $ (u) over a chosen set 
of directions fi not containing u . For equispaced linear arrays the cri- 
terion (30) leads to the Chebyshev design procedure first described by 
Dolph. However in general it is a much harder mathematical problem to 
maximize T(a) than to maximize G(a), and for that reason we shall not 
employ T{a) here. 

5. RESTRICTIONS ON SUPERDIRECTIVE ARRAYS 

Even when a set of excitation coefficients maximizing the gain function 
G(a) has been found, there is no assurance that this set will be a satis- 
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factory one on which to base the construction of a physical array. There 
is a further restriction on the solution: the array must be excited in such a 
way that when it is constructed it is likely to have a pattern $(u) which 
differs from the nominal pattern $o(u) by an acceptably small amount. 
From Section 3 the expected power pattern is of the form 



* u (u) + AV(u) £ I a, I'- 



(31) 



where A~ includes the effects of both excitation and position errors. The 
background power level relative to the main lobe of the nominal pattern 
is then just 



A'V(u) E I a k 

A— 1 

#0(110) 



a¥(u) 



K(a) 



where 



2 I a* 
K(a) = 



X) a* exp (ifJr k -Uo) 



(32) 



(33) 



K(a) is a function measuring the susceptibility of the pattern to random 
errors in the excitations and positions of the elements. Since in practice 
A J is never zero, an array with too large a value of K(a) will be unac- 
ceptable. 

Although the function K(a) has been introduced as a result of statisti- 
cal considerations, it can also be interpreted in terms of the efficiency of 
the array as an energy radiator. If we imagine the elements to have a 
certain ohmic resistance, and the excitation coefficients to correspond 
to the clement currents, then 2 [ a* | 2 is a measure of the power which is 
lost in the form of heat, and K(a) is proportional to the ratio of dis- 
sipated power to power density in the direction Uo . Thus a large value 
of K(a) corresponds to large circulating currents in the array, and to high 
ohmic losses for a given rate of radiation of power in the direction of the 
main beam. 

As a simple example of arrays which can have arbitrarily high gain at 
the expense of large values of K(a), let us consider an end-fire array of 
length L pointing in the direction u H (say the direction of the z-axis). 
The array will have n + 1 elements situated at 2 = 0, L/n, 2L/n, • • • , L. 
Let the expected excitation coefficient of the radiator at z = kL/n be 



a k = (-)'C, fc exp (ikfiL/n) 



(34) 
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where the C„, k 's are binomial coefficients. Then from (9) the norm of 
the nominal array factor is 

|/ ( U ) | 2 = | 1 - exp [i(iio-u + l)/3L/n] | 2 " 

= 2 2 " sin 2 " [(1 + cos 6)0L/2n] (35) 

= 2 2 " sin 2 " [(0L cos 2 V 2 B)/n] 

where 8 is the angle between the direction u and the direction of the 
array. If L is fixed and n is large, then approximately 

| / (u) | 2 « (2/3L/n) 2n cos 4 '* y 2 d. (36) 

Taking n large enough, one obtains an arbitrarily sharp beam and an ar- 
bitrarily high gain. On the other hand, 

f, fl ,i yv 2 - (2n)!__^_ (37) 

the last approximation being valid for large n. Hence for L fixed and n 
large, we have from (33), 

K(a) » (n//3L) 2n (n7r)- ]/2 (38) 

If A is a typical figure like 0.01, it is clear that this end-fire array is totally 
useless, in spite of its high theoretical gain, if the spacing L/n is much 
less than 1/0 = X/2x. 

The array just considered exhibits a characteristic feature of super- 
directive arrays. All such arrays depend for their narrow beamwidths on 
heavy cancellation between the fields of closely spaced radiators. In prac- 
tice the radiators can be adjusted with only finite precision, and the 
random errors which contribute to the background power level will not 
cancel out, 'on the average, as do the fields of the nominal elements in all 
directions except that of the main beam. Hence the background power 
level will completely swamp the main beam unless the array is designed 
with extraordinarily high precision. 

It turns out that merely maximizing the gain function G(a) of an array 
with interelement spacings of less than about a quarter wavelength is 
likely to lead to a superdirective design with an unacceptably large 
value of K (a). To get useful results, therefore, we should maximize G(a) 
subject to the auxiliary condition that K{a) is not to exceed a preas- 
signed value. A method for doing this is given in the next section. 

We note that it would appear possible on paper to design an array in a 
limited volume with high G(a) and low K(a) by using a large number 
of closely packed elements. For example, suppose we have an array with 
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a high value of G(a). Now suppose each element to be divided into m 
very nearby elements, each with l/m times the original excitation. The 
nominal field pattern, which determines the theoretical gain, is essentially 
unchanged, while K(a) is divided by m. The catch is that a really super- 
directive array with a reasonable value of A' would require a colossal 
number of elements to reduce K{a) to an acceptable value. Furthermore 
our statistical arguments have been based on the assumption that the 
excitation coefficients may be regarded as independent random variables 
when the antenna is built. If the elements are packed close together it 
seems unlikely that the excitations remain independent; then the use of 
K(a) to determine the precision which must be maintained is no longer 
justified. 

In dealing with superdirective arrays it should always be remembered 
that a superdirective antenna is a high-Q, small-bandwidth device, inas- 
much as the amount of reactive energy in the near field of the antenna 
is very large compared to the energy radiated per cycle. The resulting 
stringent physical limitations on superdirectivity have been exhibited 
by Chu. 9 We shall not, however, discuss questions of bandwidth here. 

6. OPTIMUM DESIGNS TAKING ACCOUNT OF FINITE PRECISION 

The procedure for maximizing the gain function 6(a) of a definite 
array, while requiring K(a) not to exceed a specified value, depends upon 
certain theorems which will now be stated. The proofs are given in the 
appendix. 

Theorem I . Let y. ^ and let a\p] denote the set of excitation coefficients 
satisfying the system of linear equations 

S hkjtijln] + ftOklii] = exp (— 0iv Uo) (39) 

y=i 

/«• = 1,2, • • • , n. Of all possible choices of a satisfying 

K(a) ^ K(aM) (40) 

the maximum value of G(a) is obtained when a = a\n]. 

The parameter ^ is essentially a Lagrangian undetermined multiplier, 
such as is commonly used in the calculus of variations and in the deter- 
mination of maxima or minima subject to a constraint. When designing 
a directive array we may select a reasonable value of K(a), say Ka , 
with an eye to the precision with which we expect to construct the array. 



9 L. J. Chu, J. App. Phys., 19, pp. 1163-1175, 194S. 
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What we should then like to do is choose the value of /x for which 

K(a{ji]) = K (41) 

and solve (39) with this value of p. Unfortunately we cannot determine n 
simply or directly from the condition (41), so we may have to make sev- 
eral trials with different values of n. If our first guess does not yield a 
value of K sufficiently close to K , the direction to proceed is indicated 

by 

Theorem II. Both K(a[n\) and G(a[n]) are monotone nonincreasing 
functions of y.; that is, dK(a\n])/dn g and dG(a\n])/dp ^ 0. 

Fairly simple expressions may be obtained (see equations (A4), (A5), 
and (A20) through (A24) of the appendix) for a{/i], K(a[n]), and G(a\ji]) 
in terms of the eigenvalues and eigenvectors of the Hermit ian matrix 
(hjk), and the parameter ju. To make a thorough study of any particular 
array configuration, one might well start by computing the eigenvalues 
and eigenvectors. Then it would be easy to plot K and G against n over 
any desired range. 

A restriction on the possible values of K(a) is given by 

Theorem III. For any array with n elements, 

K(a) ^ 1/n (42) 

and the coefficients which yield the value 1/n arc 

a, = exp (-ipik-VLo), fc = 1,2, • • • , n (43) 

up to a constant proportionality factor. 

The choice of coefficients (43) means that all elements have equal 
amplitudes and the phases are chosen to make their contributions add in 
phase in the direction u . This may be called the normal excitation; it is 
often adopted in practice as a means of obtaining a beam in the direction 
u . The patterns so obtained are most insensitive to random errors in 
the excitation coefficients, although frequently at the expense of rather 
disagreeable side lobes. Since, as pointed out in Section 5, K(a) also 
measures the efficiency of the array as an energy radiator, Theorem III 
shows that the normal excitation produces the highest power flow in the 
direction Uo for a given rate of heating the elements. 

If no auxiliary condition is imposed on the value of K(a), we have 
Theorem IV. The coefficients which maximize the gain function G(a) 
absolutely are obtained by putting n = in equations (39), so that a = a[0], 
where 

Z httaM = exp (-tfr fc -uo) (44) 

k = 1,2, •■• , n. 
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Equations (44) are equivalent to those given by Bloeh, Medhurst, and 
Pool 10 as a result of a differently phrased argument. They could also be 
obtained by following up Uzkov's remark" that the problem of max- 
imizing the directivity of an arbitrary array amounts to the problem of 
choosing an orthogonal basis for a complex linear vector space generated 
by the patterns of the individual radiators. In principle both these papers 
show how to determine the maximum gain in the general case when each 
element of the array has a different directivity pattern; but they say 
nothing about precision requirements. 

It is not invariably true that maximum gain is incompatible with 
minimum sensitivity to random errors. For example, the normal excita- 
tion of Theorem III, which minimizes K(a), simultaneously maximizes 
the conventional gain of any array of isotropic elements in which the 
distance from every element to every other element is an integral number 
of half wavelengths. Under this condition (29) shows that hjk = 5y* ; 
and the values of G and K obtained from the solution of (39) are inde- 
pendent of n. An example is a linear array with half-wavelength spacing 
between adjacent elements; another is the three-dimensional configura- 
tion of four elements at the vertices of a tetrahedron of edge A/2. If the 
elements of these arrays are not isotropic, then the normal excitation 
does not give quite the maximum gain ; but the difference will be small if 
the beamwidth of the array factor is narrow compared to the bcamwidth 
of the element pattern. For a general array with interelement spacings 
much less than a half wavelength, however, there may be a great deal of 
difference between the excitation and pattern for minimum sensitivity 
(n = oo ) and the excitation and pattern for maximum gain (/x = 0). 
Illustrative examples are worked out in Section 8. 

7. SYMMETRIC ARRAYS 

To find optimum excitation coefficients for an n-element array by the 
method of the preceding section requires the solution of a system of n 
simultaneous linear equations, and this can be laborious if n is large. 
Fortunately several types of symmetry are commonly found in antenna 
arrays. Some of these symmetries force very simple relationships to hold 
among certain of the optimum excitation coefficients, and can therefore 
be used to obtain an immediate reduction of the order of the system of 
equations (39). 

By a symmetry S of an antenna configuration we shall mean any com- 
bination of translations, rotations, and reflections which leaves the con- 

10 Reference 6, p. 304, equation (6b). 

11 Reference 5, p. 37. 
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figuration invariant. Let Sk denote the number of the element of the 
array into which the feth element is carried by the symmetry S. Some 
simple symmetries and the corresponding relationships which can exist 
between a k and a sk are listed in Table I. These relationships hold for the 
solutions a k [n] of equations (39) if the origin of the coordinate system is 
symmetrically placed,* i.e., if the origin is invariant under the symmetry 
operation S. In the last four entries of the table we assume that the pat- 
tern s 2 (u) of a single element and the weight function w(u) appearing in 
the definition of the generalized gain depend only on the angle between 

u and Un . 

To illustrate the method of proving these relationships consider the 
second one. Reflection in a plane normal to u carries a point r into 

Si = r - 2(r-u„)uo (45) 

and for any points x and y 

Uo-Sx = -Uo-x (46) 

Sx-Sy = x-y (47) 

From (46) it follows that 

s\-Su)iv(-S\i) = s\\i)w(u) (48) 

since by assumption s 2 (u) and w(u) depend only on the angle between 
u and Uo , that is, on cos _1 (u-Uo). From the last three equations and (27) 
we conclude that 



IhkM = -i- ( exp [-#(r«, - r sl! )-Sn]s\-Su)w(-Su) dSl 
4ir Js 

= -L [ exp [-ip(n -T k )-u)s\\i)w(u) dQ. 

47T Js 



(49) 



= huT 

We also have 

exp (-0r flfc -Uo) = [exp (-tj8r*-Uo)]* (50) 

The system of equations (39) can be written in the form 

X hsk.sjCisj + pa ak = exp (— «j8r s *-Uo) (51) 



fml 



* If the origin is not symmetrically placed it may be necessary to multiply 
a\n] by a suitably chosen' complex constant in order to get a set of excitation 
coefficients which satisfy the relationships of Table I. 
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or, using (49) and (50), 



2 h h j*a aj + na ak = [exp (—ipr k -Uo)]* 



or finally, taking complex conjugates, 



S hjasj* + nask* - exp (— ■ i/3r t -u ) 
j-i 



(52) 



(53) 



Comparing (39) and (53) we see that at and a 3 k* are solutions of the same 
system of linear equations. The determinant of the system is not zero, 
by equation (A7) of the appendix, and hence 

a k = a S k* (54) 

Proofs for the other cases in Table I are similar. 

In cases 3 and 5 of Table I the relationship a S k = a* reduces the system 
of equations (39) to a system of order equal to the number of different 
classes of symmetric points in the array. In cases 1, 2, and 4 the relation- 
ship a S k = dk* is more difficult to exploit. If hjk is a real matrix, which 
will be true whenever 

s 2 (-u)id(-u) = s 2 (u)w(u) (55) 

the system (39) may be rewritten as the following pair of real systems: 



X) hj Re aj 4- /z Re a* = cos (jSr* • u ) 



j-i 



2 h k j Im Oj + /x Im o fc = —sin (/3r fc -Uo) 



Now using 



Re a S j = Re a } , Im a S j = — Im oj 
the order of these systems reduces to about n/2 each. 

Table I — Relationships Among Optimum Excitation 
Coefficients in Symmetric Antenna Arrays. 



(56) 



(57) 



Case 


Symmetry 5 


Relationship 


1 
2 
3 
4 
5 


Reflection in a point 

Reflection in a plane normal to tio 

Reflection in a plane parallel to Uo 

180° rotation about an axis perpendicular to Uo 

Rotation about an axis parallel to Uo 


ask = a k * 
ask = a k * 
a S k = a k 
a S k = a** 
a S k = ak 
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It may be noted that the relationships listed in Table I still apply if 
instead of G(a) the function to be maximized is the main-beam to side- 
lobe ratio mentioned at the end of Section 4, namely 

T(a) = -JS^ m 

max $ (u) v ' 

a 

subject to the auxiliary condition on K(a). In this case we assume that 
the set S2 of directions in which radiation is to be kept small has rota- 
tional symmetry about the direction u , that is, if Ui is in fi any other 
u for which u-Uo = Ui-Uo is also in fi. 

8. NUMERICAL EXAMPLES 

We shall apply the preceding theory to two array configurations which 
illustrate some points of interest without requiring unnecessarily heavy 
computation. For simplicity throughout this section we consider only 
isotropic elements and adopt the conventional definition of gain, so that 
the matrix elements h jk are of the form given by equation (29), namely 
(sin (3rj k )/t3r jk . 

The first example consists of four elements situated at the vertices of a 
tetrahedron of edge length I (Fig. 1). The unit vector Uo is taken to be 
along the altitude of the tetrahedron which passes through the element 
a 4 and is perpendicular to the plane containing a 1} a 2 , and a 3 . This array 
is invariant under rotations of 120° about an axis containing Uo , and so 
by Case 5 of Table I we have 

ai = a 2 = a 3 (59) 

It is easily shown that the system of equations (39) for the optimum 
coefficients reduces to the two equations : 



[ 



. 2 sin /3Zl . sin 02 . 

*™& a x + (1 + m)«4 = exp [-i{%) m m 
pi 

After solving (60), the functions G(a) and K(a) may be computed from 
their definitions (26) and (33). 

The optimum gain G is plotted against the sensitivity function K in 
Fig. 1 for / = X/2, A/4, X/8, and X/16. As a matter of interest a few 
values of the parameter n are shown on the curves. The curve for I = X/2 
consists of the single point K = }i, G - 4. This is an array with half- 
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Fig. 1 — Gain G versus sensitivity K for optimum tetrahedral arrays of various 
cdnc lengths I. Some values of the parameter ju are shown on the curves. 



wavelength interelement spacing, for which, as noted at the end of Sec- 
tion (>, minimum sensitivity and maximum gain are obtained simul- 
taneously. The lower ends of the curves (/i = x) all represent normal 
excitation and minimum sensitivity; in accord with Theorem III the 
limiting value of K is }■£. The corresponding gain decreases as the size of 
the array decreases. As the condition on K is relaxed by decreasing the 
parameter n, the gain increases to a value very near 4 in all cases, and 
the sensitivity increases to larger and larger values for the smaller arrays. 
The maximum values of G, which occur for n = 0, and the corresponding 
values of K are listed in Table II. 

As a second example we consider a four-element equispaced linear 
end-fire array with interelement spacing 7 (Fig. 2). This array is sym- 
metric with respect to its center point, and so by Case 1 of Table I, 



ci\ = Qi 



(61) 



Equations (39) break up into two pairs of real equations of the form 
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Table II — Maximum Gain G and Corresponding Sensitivity K 

for Four-Element Tetrahedral and Linear Arrays with 

Interelement Spacings //A Wavelengths. 



/A 


Tetrahedral 


Linear 


G 


A 


G 


A 


y$ 

He 



4.000 
3.960 
3.990 
3.997 

4.000 


0.2500 
0.5405 
1.901 

7.371 

00 


4.000 
12.77 
15.21 
15.80 
15.95 
16.00 


0.2500 
2.065 

1.07 X 10* 
6.6 X 10 3 
4.2 X 10 B 

OO 



(56), namely: 

sin 3/3/ 



[ 



1 + M + 



3/3/ 
sin 2/3/ 



Re a x + 



sin 2/3/ 



[ sin gj sin 2/3/ 1 

L & m J 



Re ai -j- 




and 



1 + M - 



sin 3/3/ 
3/3/ 

sin 2/3/ 



Im di + 



[ sin /3/ _ sin 2/3/ 1 
L /3/ 2/3/ J 



Re (h = cos 3/3Z/2 
Re eta = cos /3//2 

[sin /3/ sin 2/3/1 T • ooi/o 

-/-^rJ Ima2= - sinW2 

[ sin /3Z~| 



(62) 



(63) 



Im ai + 1 + /x — 



Im a 2 = — sin /3//2 



/3/ 2/3/ J "' ' L" ' *" ^ 

Plots of G against K, as computed from the above equations, are 
shown in Fig. 2 for arrays with I = X/4, A/8, X/16, and X/32. Note that 
the scale of Fig. 2 is different from the scale of Fig. 1. The point for 
/ = A/2, which again falls at K = }i, G = 4, is not shown because it 
happens to coincide with the lower end of the curve for / = A/4. For 
spacings less than A/2 the curves rise to values of gain greater than 4, 
the limiting value of gain for infinitesimal spacing being 16. However, 
for spacings less than A/4 the sensitivity becomes very large for the 
higher gains, so that the upper ends of the curves are far off the hori- 
zontal scale in Fig. 2. The theoretical maximum gains and corresponding 
sensitivities are listed in Table II. 

A considerable difference is evident in the behavior of small tetrahedral 
and small linear arrays, where by "small" we mean interelement spacings 
appreciably less than a half wavelength.* Sensitivity considerations 

* It is probablynot meaningful to make a direct comparison bet\veen tetrahedral 
and linear arrays of the same spacing I, on account of their very different propor- 
tions. For example, the tetrahedral array will fit into a sphere of diameter 1.225J, 
while the linear array requires a sphere of diameter 31. 
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4 5 6 

SENSITIVITY K 



Fig. 2 — Gain G versus sensitivity K for optimum four-element linear end-fire 
arrays with various interelement spacings /. Some values of the parameter n arc 
shown on the curves. 



being ignored for the moment, the small linear array can have much 
higher gains in the end-fire direction than the tetrahedral array can have 
in any direction. On the other hand, by a suitable choice of excitation the 
tetrahedral array can exhibit very nearly the same gain in an arbitrary 
direction, ^-hile in directions different from end-fire the small linear array 
cannot produce nearly so high a gain. These observations serve to intro- 
duce a more general result, as follows: 

Let T(u) be the maximum (conventional) gain which can be achieved 
in the direction u with an array of n isotropic elements, located at ar- 
bitrary preassigned positions in space. The value of T(u) and the excita- 
tion necessary to obtain it will of course depend upon u. The following 
theorem is proved in the appendix: 

Theorem V. The average value over all directions in space of the maxi- 
mum gain T(u) of an array of isotropic elements is equal to the number of 
elements in the array; that is, 



~ f T(u) dti = n 
4ir Js 



(64) 
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From this theorem it is obvious that the tetrahedral array, which 
comes as close to spherical symmetry as is possible with four elements, 
will have a maximum gain of approximately 4 in every direction, while 
the small linear array, which has a maximum gain of nearly 16 in the 
end-fire direction, must have a maximum gain much smaller than 4 in 
most other directions. In the limiting case of an infinitesimally short 
linear array of n elements, it may be shown by a simple extension of 
Uzkov's derivation 12 that the maximum gain which can be obtained in 
the direction u is 

T(u) = E (2fc + l)P fc 2 (cos 6) (65) 

where P*(cos 6) is a Legendre polynomial and is the angle between u 
and the direction of the array. r(u) is equal to n when = 0, but its 
value decreases rapidly as varies away from zero. It is easy to verify 
Theorem V in this special case, since we have 

i- ( r(u) tffl = i E (2fc + 1) ( Pk (cos 0) sin M 

»-l ^ ' 

= I E (2k + l)[2/(2fc + 1)] = n 

An infinitesimal array is not in itself of practical importance, since for 
gains greater than unity the excitation coefficients and the sensitivity 
are infinite; but the above example does illustrate that if an array can 
have an abnormally high gain in certain directions, there must be other 
directions in which its maximum gain is abnormally low. This fact may 
have applications to the design of arrays whose physical orientation is 
fixed, while the main beam is caused to scan by electrically varying the 
excitations of the various elements. 

In conclusion we should like to point out that the design procedure 
given in this paper applies to any array in which the number of elements 
and their configuration are specified in advance. In many practical cases 
the array to be used will be so specified, and our method then yields the 
optimum design. One could, of course, ask a more general question: 
Given n elements and a sphere of fixed diameter, how should the ele- 
ments be placed inside the sphere and how should they be excited to 
maximize the gain in a specified direction, while holding the function 
K(a), which measures the sensitivity of the pattern to random varia- 
tions, below a preassigned value? The restriction to a definite number of 

12 Reference 5, pp. 37-38. 
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elements is essential, since as pointed out in Section V we could produce 
on paper an arbitrarily high gain with an arbitrarily small K(a) by pack- 
ing colossal numbers of elements into the sphere. It may well be true 
(though a proof seems to be lacking) that if there is no restriction on 
pattern sensitivity, the maximum gain possible with n elements is ob- 
tained by arranging them in an end-fire array of infinitesimal length; 
but the excitation coefficients and K(a) are then infinite. It is not self- 
evident that the end-fire arrangement will be optimum if K(a) is required 
to have a reasonable, finite value. The very interesting problem of deter- 
mining optimum array configurations for finite values of K(a) appears, 
however, to be considerably more difficult than the problem of determin- 
ing optimum excitations for arrays of fixed configuration. 

Appendix 

proofs of optimization theorems 

Let // denote the Hermit ian matrix (/&/*), a the n- vector whose com- 
ponents are the excitation coefficients, and e the n-vector whose com- 
ponents are 

o, = exp (-t/3r*-Uo), k = 1, 2, • • •, n (Al) 

The product Ha of the matrix 7/ and the vector a is the vector whose 
jth component is 

n 

(Ha) j = S h jk a k , j = 1, 2, • • ■, n (A2) 

and the inner product (x, y) of two vectors is defined as 

(as, y) = Z x k y k * (A3) 

The functions G(a) and K(a) defined by (20) and (33) may accordingly 
be written 

G(a) = L^j 2 (A4) 

{Ha, a) 

K(a) = T ^- (A5) 

I (a, e) p 

Theorem I. Consider the system of equations 

Ha + na = e, (A6) 
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where /x is a scalar parameter. If m ^ 0, then for any nonzero a whatever, 

(Ha, a) + y.(a, a) = ~ f \ /o(u) | V(u)w(u) dQ + n(a, a) > (A7) 

so that the Hermitian form (A7) is positive definite and 

det (hts + n8 kj ) ?* 

Then the system of equations (A6) does have a solution a[p]. 
From the complex conjugates of (A6) and (A7) we find that 

(aM, e) ** (A8) 

so that a[n] is not orthogonal to e. Hence an arbitrary vector a may be 
written in the form 

a = aa[n] + b (A9) 

where a is a complex scalar and 6 is a vector orthogonal to e, so that 

(b, e) = (A10) 

Note that if a were zero, we should have 

/o(uo) - (b, e) = (All) 

corresponding to no radiation in the direction u ; so we may assume 
a^O. Since G(a) and K(a) are quotients of Hermitian forms in the 
components of a, their values do not change when a is multiplied by a 
constant. Hence it suffices to consider excitations of the form 

a = afo] + b (A12) 

where b satisfies (A10). Such excitations leave (a, e) constant, so we need 
only study the behavior of (a, a) and (Ha, a) as functions of b. 
To satisfy the condition 

K(a) ^ K(a[n]) (A13) 

it is necessary that 

(a[fx), b) + (b, a\jx\) + (b, b) ^ (A14) 

Then we have 

(Ha, a) = (Habx], afc]) + (Ha[»], b) + (Hb, a[/x]) + (Hb, b) 
= (Ha\A, a\A) + (ff4»]> b) + (b, Ha\p]) + (Hb, b) 

= (Ha[n], a\n]) + (e, b) - n(a\jx], b) 

(A15) 

+ (b, e) - n(b, aM) + (#&, &) 
^ (ffaM, oM) + m(6, h) + (fffe, 6) 
^ (Ha\fi], oM) 
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where the second step follows because H is Hermitian, the third step on 
account of (A6), and the fourth step on account of (A10) and (A14). The 
equality signs in the last two steps hold if and only if b = 0. Hence the 
only possible solution to the problem of maximizing G(a) under the con- 
dition K(a) ^ K(a\p]) is given by a = aa{n], where a is a complex scalar. 
This proves Theorem I. 

Theorem II. Since H is a positive definite Hermitian matrix, it has n 
positive real eigenvalues Xi , X 2 , ■ • • , X„ , and n linearly independent 
eigenvectors v w , v {2) , • • • , v (n) , which satisfy 

Hv {r) = \ r v {r) , r = 1, 2, • • • , n (A16) 

The eigenvectors may be taken as orthonormal, so that 

(v {r \ v U) ) = 5 f( (A17) 

We may expand the right-hand side of equations (A6) in terms of the 
vectors v lr) ; thus 

n 

e = Z E r v (T) (A18) 

r=l 

where 

E r = (e, i/ r) ), r = 1, 2, • • • , n (A19) 
Assuming a similar expansion of the solution of (A6), write 

n 

alu] = Z c r v lr) (A20) 

r-l 

Substituting (A18) and (A20) into (A6) and using (A16) and (A17), we 
find without difficulty, 

c r = E r /(\ r + M ) (A21) 

We can now write down expressions for the Hermitian forms occurring 
in G(a) and K(a). Using (A18), (A20), and (A21), and the orthogonality 
condition (A17), we obtain 

n 

(flM, e) = E|£r |7(Xr + n) (A22) 

(a[/i], a[/i]) = E \E T |7(Xr + nf (A23) 

r=l 

(HahA, oM) = t, Xr | E r | 2 /(X r + M ) 2 (A24) 
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Note that all three expressions arc positive real-valued, and that 

d(a[n], e)/dn = -(«ImL afci]) (A25) 

From the preceding expressions and the definition (A5) of K(a), we 
have 

dKM) = -i-, ["(«&,], e) < g W»'frP + 2(aU a[ M ]) 2 ] (A26) 

dn (a[/x], e) J L rf M J 

The first factor on the right is positive, and the terms in square brackets 
become 

.atJ^LSJAL + ariJALTio (A27) 

r=l (X r + M) >-l (X. + M) 3 Lr=l (X r + M) 2 J 

by Schwarz's inequality. Hence K(a\ji]) is a nonincreasing function of p. 
Now let M2 > Mi and let ifi , K 2 , Gi , G 2 denote K(a[m}), K(a[n 2 ]), 
G(a\pn]), G(a[n 2 ]). By Theorem I 

G 2 = max* ( „) S jc, G{a) S max K( „)g Kl (7(a) = & (A28) 

since we have just proved K 2 ^ K\ . This completes the proof of Theorem 

II. 

Theorem III. To find the minimum possible value of K(a), let 

a = e + b (A29) 

where 

(b, c) = (A30) 

The coefficient of e may be taken as unity, since if it were zero, K(a) 
would be infinite. When a is given by (A29), (a, e) has the constant 
value 

(a, e) = n (A31) 

and making use of (A30), 

(a, a) = n + (6, 6) ^ n (A32) 

Referring back to (A5), Ave see that X(a) achieves the minimum value 
\/n when a = c; this is Theorem III. 

Theorem IV. To maximize G(a) when there is no restriction on K(a), 
set n = in the system of equations (A6). We do not require (A13) or 
(A14), but from (A 15) we still get 

(Ha, a) = (Ha[0], a[0]) + (e, b) + (b, e) + (Hb, b) 

= (Ha[Q), a[0]) + (Hb, b) (A33) 

^ (Ha[0), a[0]) 
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which shows that 6(a) is maximum when a = o[0], and thus proves 
Theorem IV. 

Theorem V. Setting p =* in (A23) and (A24) and using (A19), we 
obtain from (A4) an expression for the maximum gain in any direction 
u, namely 

T(u) = ± xr 1 \Er\ 2 = ttt \rW r) *e,V r) (A34) 

r=l r=l s=l J=l 

The dependence on u appears in the components of e, where 

e k = exp (-#r*-u), h = 1, 2, ■ • • , n (A35) 

Averaging over all directions in space gives 



— / c a ei 
-iir Js 



da = ^^' - h st (A36) 



the identification with the matrix element h, t being made from equation 
(29) on the assumption that the elements of the array are isotropic and 
the conventional definition of gain is used. It follows now from (A16) 
and (A17) that 



l/'rCuW^ZEZxrV^*/^ 

47T Js r=l *»1 (=1 

n n 

= 22 Xr _I 5 rs X, = n 

r = l .■ 1 

which proves Theorem V. 



(r) 

(A37) 



